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A Single Scalar Field Model of Dark Energy with Equation of State Crossing — 1 
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In this paper we study the possibility of building models of dark energy with equation of state 
across — 1 and propose explicitly a model with a single scalar field which gives rise to an equation of 
state larger than —1 in the past and less than —1 at the present time, consistent with the current 
observations. 


PACS number: 98.80.Cq 

Since two groups 011 independently showed the ac¬ 
celerating expansion of our universe based on the Type la 
Supernova (SN) observations of the redshift-distance re¬ 
lations in 1998, numerous observations and experiments 
consistently indicate that our universe is dominated by 
a component with a negative pressure, dubbed dark en¬ 
ergy in the literature. Especially the first year Wilkinson 
Microwave Anisotropy Probe (WMAP) data 5 strongly 
support for a concordance model in which about 73% of 
the energy of the universe is indwelt in the dark energy 
ingredient. One essential to understand the nature of 
the dark energy is to detect the evolution of its equa¬ 
tion of state w = p/p, where p and p are the pressure 
and the energy density of the dark energy respectively. 
Recent analysis on the data from the Supernova, cosmic 
microwave background (CMB) and large scale structure 
(LSS) show that the cosmological constant fits well to 
the data 0, however current data also mildly favor an 
evolving dark energy with an equation of state w larger 
than —1 in the past and less than —1 today 111 and 
evolving across —1 in the intermediate redshift. If such 
a result holds on with the accumulation of observational 
data, this would be a great challenge to the current mod¬ 
els of dark energy. Firstly, the cosmological constant as a 
candidate for dark energy will be excluded and dark en¬ 
ergy must be dynamical. Secondly, the simple dynamical 
dark energy models considered vastly in the literature 
like the quintessence or the phantom can not be satisfied 
either. 

In the quintessence models msm, w is evolving 
however always in the range of — 1 < w < 1. In the 
phantom models EL one has w < — 1. For the general 
k-essence models C = p((f, X) with X = 1/2 S7 
and the covariant derivative, it was shown that it 
also fails to cross the boundary of w = —1 because of the 
problems of singularity and classical instability in the 
perturbative equations :3i. Thus it is interesting to 
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ask if theoretically there exists model of a single scalar 
field which has a w crossing —1 during the evolution of 
the universe 1 . In this paper we propose a toy model of 
this type. 

In the general k-essence models, the Lagrangian is a 
function of the field cf> and its first derivative In 

this paper we will propose a dark energy model which 
includes higher derivative operators in the Lagrangian, 
for example, 

C = £(</>, X, □</>□</>, V^V^V^VV) , (1) 

where □ = is the d’Alembertian operator. The 

higher derivative terms in the effective Lagrangian above 
may be derived from fundamental theories. In fact it 
has been shown in the literature that this type of oper¬ 
ators does appear as the quantum corrections or due to 
the non-local physics in the string theory mum. 
With the higher derivative terms to the Einstein grav¬ 
ity, the theory is shown to be renormalizable [2^1 which 
has attracted many attentions. Recently, higher deriva¬ 
tive operators have been considered to stabilize the linear 
fluctuations in the scenario of “ghost condensation” [28} ■ 
In addition the canonical form for the higher derivative 
theory has been put forward by Ostrogradski about one 
and a half century ago f2fj. In short, it is interesting 
and worthwhile to study the implications of models with 
higher derivatives in cosmology. 

In this paper we will show this type of theories provides 
a class of models of dark energy where the equation of 
state can get across —1 during evolution, consistent with 
the current observations. We consider the simple model 
with the Lagrangian given by 2 

+ ( 2 ) 


1 The dark energy models with the property of crossing the cos¬ 
mological constant boundary are called quintom models E3 , and 
there are a lot of interests in the phenomenological stud y of this 
class of m ode ls recently in the literature Il3[ IlM Ha. Il7l Il8. 19, 

mil], mm. 

2 We take the convention for the sign of the metric with 
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where M is a constant with mass dimension and c is a 
dimensionless constant. One can see that this model dif¬ 
fers from the phantom model El due to the dimension-6 
operator. Our analysis below will show that this model 
gives rise to an equation of state which transit from 
w > -1 to w < — 1 as the redshift decreases, further¬ 
more it is equivalent to the two scalar fields model of the 
quintom dark energy studied in Ref. EL EL Ill- 

Given the Lagrangian © we have the energy- 
momentum tensor: 


= [iv^W0 + ^n0D0 + ^V^V,(n^ 


+V\gT - W0V^ - —V^V^(n^) 


( 3 ) 


and the equation of motion: 


c 2 dV 

-□ 0 -_nV + - = o. 


0) 


The energy-momentum tensor © can be rewritten as 


T = -[VpV’V^ — VpxV^x + 2V(ip — x) + 

M 2 

—xV" - + V^xV'x , ( 5 ) 

c 


where x and 4* are defined by 

c 


X M 2 D 0 ’ 

0 = 0 + X • 


( 6 ) 

( 7 ) 


It is not difficult to see that the energy-momentum tensor 
© can be derived from the following Lagrangian 


1 1 M 2 

C = --V M 0W0+-V^xWx—nV’-xK^X 2 , (8) 

with ip and \ being two independent fields. The varia¬ 
tions of the Lagrangian in © with respect to the fields 
ip and x respectively give rise to the following equations 
of motions, 


□ 

-e- 

11 

0 

( 9 ) 

M 2 


□X H-X - V = 0 , 

C 

(10) 


where the prime denotes the derivative with respect to 

0 - X- 

When we take ip = <p + x in (8) and (9), we recover 
the equation of motion Q of the single field model and 
the transformation equation ©. So, one can see from 
the derivations above that our single field model © is 
equivalent to the one with two fields © where one is 
the normal (quintessence-like) scalar field x and another 
is the ghost (phantom-like) field ip. In general, the nor¬ 
mal mode x gives rise to the equation of state w\ > — 1 


and the ghost mode has w 2 < —1, consequently the en¬ 
ergy density in x would decrease and that in ip would 
increase with the expansion of the universe. Assuming % 
dominates the dark energy sector in the early time, the 
equation of state for the system will be larger than — 1 
during that period of time. When the ghost mode be¬ 
comes dominant over the normal one (and will continue 
to be dominant among all the components in the uni¬ 
verse), the equation of state of the system will become 
< — 1 and cross the cosmological constant boundary at a 
intermediate redshift. Quantitatively the moment when 
the equation of state crosses —1 depends on the values of 
the model parameters. 

Now we study the perturbations of our model. As is 
well known, in the usual quintessence and phantom mod¬ 
els, the field and metric perturbations are stable on the 
small length scales. For both the quintessence and the 
phantom models the sound speed as indicated in Ref. 
|30j is given by 


2 __ dp/dX _ 

Cs dp/dx 


( 11 ) 


where X = (1/2 )V^cpX^cp. This can be seen clearly from 
the dispersion relations. For the phantom field, the fre¬ 
quency u relates to the momentum k in the same way as 
in the case of the quintessence, w 2 = k 2 + .... And the 
sound speed is often defined as uj 2 = c 2 k 2 + .... Since 
we have demonstrated that our model is equivalent to a 
model with a quintessence-like and a phantom-like scalar 
fields, the general arguments for the quintessence and the 
phantom hold here, which indicate there will be no expo¬ 
nentially instabilities of the perturbations on small scales 
in our model. 

In general the potential term V(ip — x) (equivalently 
V{(p) in our single field model) should include the inter¬ 
actions between the two fields ip and y. For some specific 
choices of the potential V, however these two modes could 
decouple. As an example we consider V = (1/2 )m 2 (p 2 . 
Then the Lagrangian © can be “diagonalized” as 


£ = - ^m\(p\ - ^ m 2 2 (pl 


through the transformation 


ip 
X 

In (12) and (13) 
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Hence, one can see that our model © with V = 
(1/2 )m 2 (f) 2 is equivalent to the uncoupled system 1121) . 
The two modes fa and fa evolve independently in the 
universe. The positivity of the parameter c guarantees 
the positivity of m\ (m 2 is always positive as long as 
1 + >0). In the limit of m <C M, the masses 

of the two modes are approximately m\ — M/fac and 
m -2 — m. In fact, <j>\ and fa are the eigenfunctions of the 
d’Alembertian operator, □, with the eigenvalues —m 2 
and m 2 respectively. The solution to the equation of 
motion 0 4> is decomposed by these eigenfunctions as 

0 = (1 + + fa )- 

As we mensioned above, the perturbations in the phan¬ 
tom component fa present no unphysical instabilities in 
the classical level. In fact, the spatial fluctuation of 
fa with wavenumber k > m 2 is stable, the solution to 
the equation of motion is oscillatory with time. The 
time-rising behavior only presents on very large scales, 

L > 1/to 2 . If the phantom mass m 2 of the phantom field 
is smaller than H , there will be no instability within the 
horizon [31|| . The rising behaviors of the super-horizon 
modes (including the background mode, which corre¬ 
sponds to k = 0) of phantom do not matter in cosmol¬ 
ogy. This is because the universe is expanding instead 
of static, generally the components in the universe are 
evolving. Furthermore, the Hubble expansion provides a 
friction force which prevents these modes from increas¬ 
ing exponentially. In another word, such an instability is 
“benign” 

In Fig0 we plot the equation of state of our model 
as a function of In a. One can see clearly w crosses —1 
during evolution. Furthermore, we made a test on the 
equivalence of the model 0 with the two-field quin- 
tom model numerically for V = (1/2)to 2 ^ 2 . If the 
background spacetime is described by the Friedmann- 
Robertson-Walker metric, ds 2 = dt 2 — a 2 (t)5ijdx 1 dx J 
and the field <j> is homogeneous, the d’Alembertian op¬ 
erator will depend only on the time derivatives, □ = 
d 2 /dt 2 + 3Hd/dt. The equation of motion 0 becomes: 

rj 2 O 

-(f)-3H(f) + m 2 (t) - ^ 2(^2 +3H—)((j) + 3 H (j)) = 0 , 

(16) 

where dot represents the derivative with respect to time. 
The energy density and the pressure are 

p = ~^ 2+ 2 ^^ +3H ^ )2 -; 

(17) 

P = ^0|(0+3I?0)-im 2 0 2 , 

(18) 

respectively. The initial conditions picked up in the nu¬ 
merically calculations are 

0(0) = <?[0i(0) + 02(0)] , 

0(0) = 001 (o) + 02(0)] , 

0(0) = 0-3#(O)(0i(O) + </ 2 (O))-to 2 </i(O)+TO202(O)], 



FIG. 1: Plots of the equation of state for our model and the 
two-field model with one quintessence and one phantom field. 
In the numerical calculations we choose the initial conditions 
given by Eq. m and have set m = 2 x 10 61 M p i, = 
10 120 M~[ 2 . The equivalent initial conditions for the two scalar 
fields are </>i(0) = 3 x 10 _2 M p i, c/>2(0) = 0.18M p i, and 0i(O) = 
0, 02(0) = 0. f Ide ~ 0.73 and h « 0.68. 


0(0) = 5 [(9i7(O) 2 -3i7(O))(0 1 (O) + 0 2 (O)) + 
3ff(O)(m 2 ^i(O)-m 2 0 2 (O))- 
(m^0i(O) - m|0 2 (0))] , 

H 2 ( 0 ) = ^[pm{Q) + P M + \m-\^m + 

0) + 7^202(0)] > 

77(0) = -47rG[p m (0) + ^p r ( 0) + 0?(O) - 0|(O)] , (19) 

where 0 in the bracket means the initial time, p m and p r 
are the energy densities of matter and radiation respec¬ 
tively, and g is defined as 

9 = (1 + ( 20 ) 

We depicte the numerical result for the comparison in 
FigHI One can see from this figure the equation of state 
from 03 and © of our single scalar field model does 
coincide with that of the two-field model 03 - 
In summary, we in this paper have proposed a single 
scalar field model of dark energy which has the property 
of crossing the cosmological constant boundary. This is 
achieved by considering higher derivative operators. And 
this class of model is free from the difficulties of the sin¬ 
gularity and the gravitational instabilities in the general 
k-essence-like models. We should point out that although 
the energy-momentum tensor of this single scalar field 
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model is equivalent to two independent scalar fields 3 , 
when we consider the interactions they may possibly 
show some different behaviors. It might be theoretically 
possible to solve the inherent problems, the quantum in¬ 
stabilities exhibited in the model of phantom [sdj and is 
worth studying further. 

Our study here has the following implications: 1) in the 
general phenomenological fitting appearing in the litera¬ 
ture of the dark energy models to the observational data 
such as SN, CMB, LSS an issue is how to consistently 
include the perturbations of the dark energy. The key to 
this issue is to build theoretical models of dark energy 
with the equation of state crossing — 1. Our model pro¬ 
posed in this paper serves as an example for this class 
of models of dark energy to study the cosmological irn- 

3 After our paper appeared, the authors of Ref fill used the 

term only and dropped the conventional kinetic term as well the 


plications numerically; 2) similar to the k-inflation our 
model can be applied for inflation. And in this case one 
would expect non-vanishing isocurvature perturbations 
even though this model includes explicitly only one infla- 
ton field. 
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potential, which cannot be identified with two scalar fields. 


[1] Riess AG et al, 1998 Astron. J. 116 1009 
astro-ph /9805201 

[2] Perlmutter S et al, 1999 Astrophys. J. 517 565 
astro-ph/9812133 

[3] Spergel D N et al, 2003 Astrophys. J. Suppl. 148 175 
astro-ph /0302209 

[4] Seljak U et al., 2005 Phys. Rev. D 71 103515 
astro-ph /0407372 

[5] Alam U, Sahni V and Starobinsky A A, 2004 JCAP 0406 
008 astro-ph/0403687 ; Huterer D and Cooray A, 2005 
Phys. Rev. D 71 023506 astro-ph/0404062 ; Gong Y, 
2005 Class. Quant. Gray. 22 2121 astro-ph/0405446 

[6] Corasaniti P S, Kunz M, Parkinson D, Copeland E 
J and Bassett B A, 2004 Phys. Rev. D 70 083006 
astro-ph/0406608 ; Hannestad S and Mortsell E, 2004 
JCAP 0409 001 astro-ph/0407259 ; Upadhye A, Ishak 
M and Steinhardt P J, 2005 Phys. Rev. D 72 063501 
astro-ph/0411803 

[7] Wetterich C, 1988 Nucl. Phys. B 302 668 

[8] Wetterich C, 1995 Astron. Astrophys. 301 321 

hep-th/9408025 


[9] 

[ 10 ] 

[ 11 ] 

[ 12 ] 

[13] 

[14] 


Ratra B and Peebles P J E, 1988 Phys. Rev. D 37 3406 
Caldwell R R, Dave R and Steinhardt P J, 1998 Phys. 
Rev. Lett. 80 1582 astro-ph/9708069 
Caldwell R R, 2002 Phys. Lett. B 545 23 

astro-ph /9908168 

Zhao G -B, Xia J -Q, Li M, Feng B and Zhang X, 2005 


Preprint astro-ph/0507482 

Vikman A, 2005 Phys. Rev. D 71 023515 

astro-ph/0407107 

Feng F, Wang X and Zhang X, 2005 Phys. Lett. B 607 35 
astro-ph/0404224 ; Feng B, Li M, Piao Y -S and Zhang 
X, 2004 Preprint astro-ph/0407432 

[15] Wei Y -H and Tian Y, 2004 Class. Quant. Grav. 21 5347 
gr-qc/0405038 

[16] Zhang X, 2004 Preprint hep-ph/0410292 

[17] Guo Z -K, Piao Y -S, Zhang X and Zhang Y -Z, 2005 
Phys. Lett. B 608 177 astro-ph/0410654 

[18] Hu W, 2005 Phys. Rev. D 71 047301 astro-ph/0410680 


[19] Xia J -Q, Feng B and Zhang X, 2005 Mod. Phys. Lett. A 
20 2409 astro-ph/0411501 

[20] Caldwell R R and Doran M, 2005 Phys. Rev. D 72 043527 
astro-ph/0501104 

[21] Li B J, Chu M -C, Cheung K C and Tang A, 2005 
Preprint astro-ph/0501367 

[22] Wei H and Cai R -G, 2005 Class. Quant. Grav. 22 3189 
hep-th/0501160 

[23] Zhang X F, Li H, Piao Y -S and Zhang X, 2005 Preprint 
astro-ph/0501652 

[24] Simon J Z, 1990 Phys. Rev. D 41 3720 

[25] Elizer A and Woodard R P, 1989 Nucl. Phys. B 325 389 

[26] Erler T G and Gross D J, Preprint hep-th/0406199 

[27] Stelle K S, 1977 Phys. Rev. D 16 953 

[28] Arkani-Hamed N, Cheng H -C, Luty M A and Muko- 
liyama S, 2004 JHEP 0405 074 hep-th/0312099 ; Arkani- 
Hamed N, Creminelli P, Mukohyama S and Zal- 
darriaga M, 2004 JCAP 0404 001 hep-th/0312100 ; 
Anisimov A and Vikman A, 2005 JCAP 0504 009 
hep-ph/0411089 ; Mukohyama S, 2005 Phys. Rev. D 71 
104019 hep-th/0502189 ; Arkani-Hamed N, Cheng H -C, 
Luty M A, Mukohyama S and Wiseman T, 2005 Preprint 
liep-ph/0507120 

[29] Ostrogradski M, 1850 Mem.. Ac. St. Petersbourg VI 4 385 

[30] Garriga J and Mukhanov V F, 1999 Phys. Lett. B 458 
219 hep-th/9904176 

[31] Hsu SDH, Jenkins A and Wise M B, 2004 Phys. Lett. 
B 597 270 astro-ph/0406043 ; Buniy R V and Hsu S D 
H, 2005 Preprint hep-th/0502203 

[32] Smilga A V, 2005 Nucl. Phys. B 706 598 

hep-th/0407231 ; 2005 Preprint hep-th/0503213 

[33] Anisimov A, Babichev E and Vikman A, 2005 JCAP 0506 
006 astro-ph/0504560 

[34] Carroll S M, Hoffman M and Trodden M, 2003 Phys. 
Rev. D 68 023509 astro-ph/0301273 ; Cline J M, Jeon 
S -Y and Moore G D, 2004 Phys. Rev. D 70 043543 
hep-ph/0311312 















